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Abstract 

The standard model prediction for a magnetic dipole moment of a neutrino is pro- 
portional to the neutrino mass and extremely small. It also generates a flavor-changing 
process, but the GIM mechanism reduces the corresponding amplitude. These proper- 
ties of a neutrino magnetic moment change drastically in a degenerate electron plasma. 
We have shown that an electron-hole excitation gives a contribution proportional to the 
electrons' Fermi momentum. Since this effect is absent in /i and r sector, the GIM can- 
cellation does not work. The magnetic moment induces a neutrino oscillation if a strong 
enough magnetic field exists in the plasma. The required magnitude of the field strength 
that affects the v e burst from a supernova is estimated to be the order of 10 8 Gauss. 
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A magnetic dipole moment of a neutrino induces interesting phenomena, such as a spin 
rotation when its travelling in a static magnetic field [0], |2j or a transition radiation when 
passing an interface between two different media ||. It may also affect the stellar cooling by 
the decay of plasmons H, which is known as a dominant cooling process in a dense star || ||. 
The standard model predicts a nonzero value for it through the processes depicted in Fig. 1, 
but the magnitude is far below the one that current experiments can detect. To leading order 
in mf/rriyy, the result is independent of the mass m>i of the charged leptons in the internal lines; 
it is instead proportional to the neutrino mass m v @, 



where hb is the Bohr magneton. There are many orders of magnitude between this prediction 
and the present empirical limits, \i v < 1CT 6 — 10 _11 /iB |§. 

The leptonic charged current coupling to the W boson generates a generation mixing for 
massive neutrinos as is the case for the quark sector. The processes in Fig. 1 can then gen- 
erate, so to say, a flavor- changing electromagnetic current. The independence of the leading 
contribution from m/, however, subjects the corresponding amplitude to the GIM suppression 
[0] : The sum of the leading contributions from all three generations cancels each other because 
of the unitarity of the leptonic CKM mixing matrix. Thus the amplitude, such as the one for 
the decay — > z/ e 7, gets a nonzero contribution from the next-to-leading effect and is further 
suppressed by a factor (mf/m^) than one estimates naively with \i v . 

The purpose of this letter is to show that these properties of a neutrino magnetic moment 
change drastically in a degenerate electron plasma. In the gas of high density and relatively low 
temperature, i.e., where the Fermi momentum p-p is much larger than the temperature T, most 
of the electrons degenerates into the Fermi sphere. Electromagnetic potential excites one of 
the electrons out of the Fermi sphere and leaves a hole in it. Subsequently the excited electron 
comes back into the sphere emitting a pair of neutrinos by exchanging a W boson with the 
hole. This polarization of a electron-hole pair turns out to add a contribution of the order of 
cGyPf to the magnetic moment (See Eqs. ([11]) and (0) below). Since this effect is intrinsic to 
degenerate electrons and absent in the [i and r sectors, the GIM mechanism no longer washes 
out this contribution. 




(1) 
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The Fermi momentum of the electron gas at the core of massive stars becomes as large as or 
larger than the electron mass m e at the later stage of their evolution. This gives a possibility 
that the resulting magnetic moment becomes so large that some observations can detect it. 
We will consider a neutrino oscillation induced by the flavor-changing magnetic moment, which 
may take place in the stars. Our calculation uses a zero temperature approximation, and the 
result should be applied to the cases of pp ^> T. In the following, we will first describe the 
result for the magnetic moment briefly and then consider the induced neutrino oscillation. 

We assume the masses that the neutrinos, v e , u^, and u T , get at the electroweak symmetry 
breaking are of Majorana type.[] The corresponding fields, which we collectively denote by v{x) 
in the left-handed two-component notation, have an expansion 

v{x) = \ E [e- iE ^ x ° +ip - 3 u(p, s) afriy) + e^>°~^ v{p, s) aj.»] (2) 

in the annihilation, a^ s (u), and the creation, at s (V), operators of the state with the momentum 
p and the helicity s; the spinors u(p, s) and v(p, s) are defined by 



u(p, s) 
v{p, s) 



E(p) — s\p 



xip, s) (3) 



Jn e X (p,s)*, (4) 



E(p) + s\p\ 



\ 2E{p) 

where E(p) = (p 2 + m^ 2 ) 1 ^ 2 , s) is the normalized eigenspinor for the helicity, 

(p- a)x(p,s) = sx(p,s), xHP, s )x(P,t) = S st , (5) 

defined with the Pauli matrices a, and e = ia 1 . The quantisation volume V in (fj) should 
be taken as the size of the plasma. Note that for states with relativistic momentum, v(x) is 
dominated by annihilation operators of s = — 1 and creation operators of +1. 

The electrons and positrons are described by e(x) and e°(x) in the two-component notation. 
e(x) has the charged current coupling to a W boson, while e°(x) does not. They have the 
same expansion as Eq. @ if one does proper replacements; in e(x) the annihilation operators 
stand for electrons, a^ s (e), and the creation operators for positrons, &^ s (e), while in e c (x) the 

annihilation operators stand for positrons, b^ s (e), and the creation operators for electrons, 

1 For the case of Dirac neutrinos, one readily gets the corrections by simply replacing at s (y) in (||) with 
bp s (v), the creation operator for the anti-neutrino state. It does not affect our discussion in this letter. 



a|. s (e); the energy in the expressions of u(p, s) and v (p, s) is understood to be calculated with 
m e . 

In terms of these fields, the one-loop induced electromagnetic vertex for the neutrinos has 
the form 

*£ eff = {^f\ V„ aC V; b A(Pi)F»(p l ,P2)MP2)A tl ( q ), (6) 

where, e and g are the electromagnetic and charged current coupling constants, V Uae is the 
element of the CKM matrix between v a (a = e, fi, r) and e, and the fields are written in their 
Fourier components of momenta q, p\, and p%, that satisfy q + p\ + P2 = 0. The structure 
function in Eq. (|^) is given from the expectation value of a T-product, 

F"(p 1 ,p 2 ) = J J dxdye- ipiX - ip2V 

x(\TW+(x)a v e(x) ^(0)^(0) - e ct (0)a^e c (0) + ...] e< (y)a x W^ (y)\ ), (7) 

where the state 

i>= n 4»i°> ( 8 ) 

\p\<p F ,s=±l 

represents the degenerate electron plasma; <r M = (1, — a). We have abbreviated the contribution 
from W's electromagnetic vertex in Eq. (0). The calculation is carried out by modifying the one 
for a non-relativistic plasma 0. The usual Feynman rule applies if one uses the propagators 
that take Pauli exclusion principle into account. They are obtained from Eq. (1) in Ref. ||. 
Since we are interested in the effect proper to a degenerate plasma, we subtract the vacuum 
contribution F£. from F M and concentrate on the remaining term F£ = F M — F£. (F^ c is defined 
by F p with zero pp). 

We first specify the form of F£ taking various conditions into account. Since we have 
assumed the plasma is isotropic and homogeneous, the temporal component F° g is a scalar 
while the spacial components F^ g are vectors. To leading order in l/m\v, only the contribution 
of Fig. 1 (a) remains and the W boson propagator can be safely contracted to a point form, 
Gfiv ~ SvA^w 2 ! - This is because the loop momentum is restricted by pf in F^ g . Thus F£ 
depends only on q. There are two structures, 1 (unit matrix) and a ■ q, for and four 
structures, q, (<? • q)q, a, and axg, for F dg . They are also constrained by the gauge invariance, 
q°F® — q ■ Fd g = 0. Finally, we recall the study on F£ in Ref. for the plasmon decay. 
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Although it was done with the four Fermi interaction, the result still applies because the W 
boson propagator is safely contracted. The four Fermi interaction is Fierz-transformed into the 
product of the neutral current of neutrinos and the V-A current of electrons. F^ g is then divided 
into two components, one from electrons' vector current and the other from their axial vector 
current |5[]. The vector component is related by F£ g ~ U^cr" /2my/ 2 to the polarization tensor 
IT^i, and, thus, turns out to consist of two independent structures. For the axial component, 
we are left with only one available structure a x q. We realize that F£, has the form 



,0 



fi q"q 



1 + 



q 



\Q\ 



+ ft W\ 2 v - {S • q)q) +/ m i(ffx? 



(9) 
(10) 



with three form factors f\, ft, and f m , which depend on the rotational invariants \q \ and 
Applying the time reversal invariance of | ), we see that all of the form factors are an even 
function of q . 

The flavor- changing magnetic dipole moment is related to f m by 

,2 



H>ab — 2 V v a cV Ub cJrcii 

where we have used the definition that the corresponding Lagrangian is written as 



(11) 



(12) 



with the magnetic field B. Under the CP transformation, / m is odd (while f\ and f t are even). 
That means fi a b changes its sign if the plasma is made of positrons, or positrons and electrons 
contribute destructively to fi a b if they coexist in a plasma. 

We evaluated I = Tr(cr : 'F ( | g )/2 to extract out / m . Keeping only the terms that are propor- 
tional to the anti-symmetric tensor e^, we found 

1 \ f dp 



V 



q k _ P k \ 0(\p + q\~ Pf)0(pf - \p\)(E(p + q)- E(p)) 
\E(p + q) E(p)J (q°) 2 -((E(p + q)-E(p)) 2 

'p k + q k p k \ 9(p F -\p + q\)(E(p + q) + E{p)) 

E(p + q) + W)) (q ) 2 -((E(P + <i) + E(p)) 2 



(13) 



2 As far as the decay of a plasmon is concerned, f\ and ft give dominant contribution over f m || (See also 
11)- 



where E{p) = {p 2 + m 2 ) 1 / 2 . The first term comes from the electron- hole excitation; the 
second represents the vacuum polarization of a electron-positron pair that is now forbidden if 
the electron has the momentum in the Fermi sphere. We are interested in the static component 
of fjLab] for q° = and \q | <C pp, we found 



/- = - J-2—2- ( 14 ) 

An 2 m-w 



Details of the calculation will be presented elsewhere JTT 



Let us turn to a neutrino oscillation induced by fiat,. The oscillation we are considering here 
is the one where a neutrino in one of the flavors (mass eigenstates) oscillates into another flavor 
in the presence of an external magnetic field B. Thus, it is conceptionally different from the 
vacuum oscillation [1^] in which a neutrino, created in an eigenstate of the weak interaction 
coupled to an electron, oscillates into another kind. We assume the neutrino energy E is 
relativistic, m v jE <C 1, in the Lorentz flame where the plasma has zero mean velocity. We also 
assume the deviation of B from the completely static and homogeneous configuration is small, 
and the energy and momentum transfer from B to the neutrino is negligible compared with E. 

In a vacuum, where a flavor-changing process is suppressed, a physically interesting process 
induced by a magnetic moment is an oscillation between two different helicities. This is the spin 
rotation discussed in Refs. |TJ, 0]- This helicity-flipping process is, however, suppressed by the 
factor (m u /E) as one can immediately see from the explicit form for the coupling, Eq. (JT^), and 
the expansion, Eq. An advantage in a degenerate plasma is that an oscillation is possible 
between two different flavors keeping the helicity intact and thus without the suppression of 
m v /E. 

For the sake of clearness of the argument, we specifically consider an oscillation between u e 
and with s = — 1. We adopt a two-state approximation; we restrict the Hilbert space with 
two states, \u e ) = i(v e )| ) an d |^) = a$ 2 _i(^)l )■> anc ^ consider an oscillation just between 
them. Based on the assumptions mentioned above, we assume the energies are the same. Then 



the momenta (whose magnitudes are yj E 2 — m 2 e and E 2 — m 2 ^) should be slightly different, 
which is taken account for by a small momentum transfer q from B. 

The Hamiltonian H, a 2 x 2 matrix in our approximation, get off-diagonal matrix elements 



3 Thc proportionality of //„ to m„ comes from this factor in our two component notation. 
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in the presence of B. It reads 

a ={i>i)< < i5 > 

where 

A= dx(v c \(— £ m )|z/ M ) ~ /i eAt cos9\B\, p-B = cos9. (16) 

The new eigenstates are a mixture of |z/ e ) and with equal weights and a different relative 
phase; their energy difference is AE = 2\A\. The half-oscillation length = n/2AE 

characterizes the oscillation; the relativistic neutrino that is initially in |z/ e ) gets a fifty percent 
probability to be detected as after a travel of this length. 

We mention a few of the characteristics of i/ 1 / 2 ) , other than the absence of the suppression 
of m v /E mentioned above. First, it does not depend on the energy or the mass difference of the 
neutrinos, while in the vacuum oscillation the oscillation length is given by 7rE/(m^ — m^). 
Secondly, it inversely proportional to cos#; the oscillation takes place most efficiently for the 
momenta nearly parallel to B, while the spin rotation is most efficient when the neutrino travels 
perpendicular to B [0]]. 

The stellar interiors are the candidates where the magnetic-moment-induced neutrino oscil- 
lation may have a physical importance. We relate Pf to the mass density p of the interior and 
the electron's number fraction Y e per baryon, and write the oscillation length as 

l/ 1 / 2 ) = 6.1 x 10 17 cm[F c p(g/cm 3 )]- 1 / 3 [| J B|(G)]- 1 [|K ee K; e cos^]- 1 (17) 

with values of p in units of gram per cubic centimeter and of \B | in Gauss. Our zero temperature 
calculation applies for the case pf 3> T, which reads [Y e p(g/ cm 3 )] 1 / 3 ^> 1.7 x 10 _8 T(K) with 
T in units of Kelvin. For the solar neutrinos, the oscillation is too slow; even taking optimistic 
values, p ~ 10 2 g/cm 3 and \B\ ~ 10 3 G, we get ~ 10 14 cmY-^HV^V*^ cos#|]-\ which 

is much larger than the solar radius ~7x 10 10 cm. 

An interesting possibility is that the oscillation may convert z/ e from a supernova to or 
v T . A massive star, whose mass is bigger than 8M , has its core consist mainly of the Fe 
elements at the final stage of its evolution. The degenerate electron plasma plays an important 
role to keep the core from collapsing. It, however, fails when the mass of the core exceeds the 
Chandrasekhar mass. The core begins to collapse and u e is copiously produced by the electron 
capture by the Fe nuclei. The collapse eventually stops when the central density reaches to 
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the nuclear density, and a shock wave is formed at the central region of the core. The shock 
wave then produces the u e burst by the neutronization of free protons when it propagates 



outward in the core |T3[. The v e flux generated by these processes in the very early stage of 
the explosion is converted to another kind if there penetrates a strong enough magnetic field 
in the core; the resulting flux becomes a composition of v e and a substantial amount of or 
v T . We estimate the required field strength for this to happen by the condition that L^ 1 ' 2 ^ is 



,3 



shorter than R COTe , the core radius. Using the values R 

core ~ lO^cm and p ~ 10^g/cm 
we find \B \ > 10 8 G ^ _1 ' 3 [|K, ee V^* e cos#|] _1 . This magnitude seems realizable if one compares 
it with the value \B | ~ 10 13 G which is possible for a supernova or a neutron star as discussed 
in Ref. 0. 
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Figure Caption 

Fig. 1: The Feynman diagrams for the magnetic dipole moment of a neutrino in the standard 
model. 
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